MPP-2004-171 



Electroweak Absorptive Parts in NRQCD Matching Conditions 



Andre H. Hoang and Christoph J. ReiBer 



Max-Planck-Institut fur Physik 
(Werner-Heisenberg-Institut) , 

Fdhringer Ring 6, 
80805 Miinchen, Germany * 



Abstract 

Electroweak corrections associated with the instability of the top quark to the next-to-next- 
to-leading logarithmic (NNLL) total top pair threshold cross section in e + e~ annihilation are 
determined. Our method is based on absorptive parts in electroweak matching conditions of the 
NRQCD operators and the optical theorem. The corrections lead to ultraviolet phase space diver- 
gences that have to be renormalized and lead to NLL mixing effects. Numerically, the corrections 
can amount to several percent and are comparable to the known NNLL QCD corrections. 
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I. INTRODUCTION 



The line-shape scan of the threshold top pair production cross section cr(e + e~ — > 7*, Z* — > 
tt) constitutes a major part of the top quark physics program at the International Linear 
Collider (ILC) project that is currently being initiated. Because in the Standard Model 
the top quark width T t « 1.5 GeV is much larger than the typical hadronization energy 
Aqcd, it is expected that the line-shape of the total cross section is a smooth function 
of the c. m. energy, and that non-perturbative effects are strongly suppressed [3, 0| • The 
determination of the top quark mass (in a threshold mass scheme |3|, |^] ) is the most important 
measurement that can be obtained from the threshold scan since an uncertainty of only 
around 100 MeV is expected 0, 0|. This prospect is quite robust, from the theoretical as 
well as from the experimental point of view, since it relies mostly on the determination of 
the c. m. energy where the cross section rises. Because the tt pair is produced predominantly 
in an S-wave state the rise of the cross section is quite rapid and easily measurable even 
in the presence of beam effects j5(. In addition it will also be possible to determine the 
strong coupling a s , the total top quark width T t and, if the Higgs boson is light, the top 
Yukawa coupling g t th- However, the latter measurements are sensitive to the form and the 
normalization of the line-shape. Since the observable cross section is a convolution of the 
theory prediction with the partly machine- dependent luminosity spectrum arising from QED 
effects [1,0], high demands are imposed on theoretical predictions and experimental analyses 
to make these measurements possible. In particular, theoretical predictions need to have a 
precision at the level of only a few percent. 

The common theoretical tool to make computations for top threshold observables is non- 
relativistic QCD (NRQCD). 1 It provides an economic and systematic treatment for the non- 
relativistic expansion and for QCD radiative corrections coming from high and low energies. 
It has also become evident that it is advantageous to use renormalization group methods [1[ 
to resum logarithms of the top quark velocity v to all orders of QCD perturbation theory in 
order to avoid large normalization uncertainties of at least 20 % that are obtained in fixed- 
order predictions p] . Concerning QCD effects at next-to-leading-logarithmic (NLL) order for 
the total cross section only the tt production current has a non-trivial running, which is fully 
known 0,0,0, 11 1. At next-to-next-to- leading- logarithmic (NNLL) order for the total cross 



section the QCD evolution of almost all required couplings is known [§, [l^, IH 3| except for 
missing subleading mixing effects in the running of the heavy quark pair production current. 
Theoretical analyses for the total cross section at NNLL order were given in 0, 0, [l^f. 
Currently the normalization uncertainties of the total cross section from QCD effects are 
estimated to be around 6% [I?}- 

While the major focus in recent work was directed on a better understanding of QCD 
corrections at NNLL order, a systematic treatment of electroweak effects at the same level 
has not yet been accomplished. Electroweak corrections are responsible for a variety of 
different physical effects. At leading order, the three basic electroweak effects are the e + e _ 
annihilation process that leads to top pair production by virtual photon and Z exchange, 
the finite top lifetime, which can be implemented by the replacement rule E — ► E + iT t [2| 



We use the term NRQCD to refer to a generic low-energy effective theory which describes nonrelativistic tt 
pairs and bound state effects and not for a theory valid only for scales m t > fi > m t v. For the presentation 
we use the conventions and notations of vNRQCD established in ||| Q , but we emphasize that our results 
are generally true. 



2 



(E = y/s — 2m t being the c. m. energy with respect to the top threshold and T t being the 
on-shell top quark width), and the luminosity spectrum mentioned above. All three effects 
have so far been treated independently. It has become the convention that only the first two 
effects are included into theoretical predictions while the luminosity spectrum is accounted 
for in the experimental simulations p . 

At the subleading level a coherent treatment of electroweak effects has not yet been 
achieved, although previous partial analyses have indicated that they can reach the level 
of a few percent [HI, 19| . It is also evident that at the subleading level an independent 



treatment of various different electroweak effects will, eventually, be impossible. Since a 
systematic treatment relies on the consistent separation of off-shell (non-resonant) and close- 
to-mass-shell (resonant) fluctuations, the concept of effective theories appears again to be 
a highly efficient tool to make progress [20]. For the ti threshold such a framework is 
already provided by the effective theory NRQCD itself, since it automatically achieves an 
expansion in the off-shellness of the top quark through the non-relativistic expansion in 
v. The NRQCD effective theory formalism can be extended to account for electroweak 
corrections. For example, including the electroweak radiative corrections to the top quark 
two-point function in the NRQCD matching conditions one obtains the additional heavy 
quark bilinear terms 2 

p p 

in the effective Lagrangian, where ijj p and x P represent Pauli spinor field operators destroying 
top and antitop quarks, respectively. For simplicity color indices are suppressed throughout 
this paper. The terms in Eq. (JQ reproduce the replacement rule E — > E + iT t . They render 
the effective theory Lagrangian non-hermitian because they describe an absorptive on-shell 
process that has been integrated out from the theory. Nevertheless, they allow for a correct 
determination of the total cross section from the forward scattering amplitude using the 
optical theorem and the unitarity of the underlying theory. In fact, this issue is analogous 
to the so-called strong phases in QCD amplitudes that are the basis for the search for CKM 
CP-violation in a number of B meson decays [2l| . 

In this paper we extend this approach and investigate the role of absorptive parts related 
to the top quark decay in electroweak loop corrections to the NRQCD matching conditions 
that contribute to the NNLL total cross section. 3 For simplicity we neglect the bottom quark 
mass and set V t b = 1, and approximate the W boson and bottom quark as stable particles. 
We demonstrate that these electroweak corrections properly account for the interference 
of the dominant double-resonant process e + e _ — ► tt — > bW + bW~ with the w 2 -suppressed 
single-resonant amplitudes for e + e~ — > bW + i — > bW + bW~ and e + e~ — > tbW~ — > bW + bW~ . 
We also show that absorptive parts that do not contribute to the bW + bW~ final state, and 
are therefore also not accounted for in line-shape measurements, can be excluded in a gauge 
invariant way. This requires to include also fields for the electrons and positrons from the 
initial state into the effective theory, which act like classic fields for QCD interactions. The 



2 Note that in the effective theory we use positive energy spinors for the antiparticles. 

3 In a complete treatment of all electroweak effects one can integrate out all electroweak effects associated 
with the massive W, Z and Higgs bosons at the scale m t . Below the scale m t gluons, photons and quarks 
remain as dynamical degrees of freedom. 
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new corrections show interesting features. They slightly modify the form of the cross section 
line-shape and lead to UV phase space divergences that are directly related to the fact that 
the top quarks are unstable. The divergences lead to anomalous dimensions of (e + e - )(e + e~) 
operators that also contribute to the absorptive part of the forward scattering amplitude. 
In the total cross section these mixing effects contribute at NLL order and represent a novel 
NLL effect that has remained unnoticed in previous work. Numerically, the size of the new 
corrections ranges up to 5% and partly compensates for the large QCD corrections found 
recently in [l3| . 

The program of this paper is as follows. In Sec.|n]we discuss the power counting for the 
total cross section with respect to electroweak effects associated with the top quark decay. 
We introduce the (e + e~)(t t) effective theory operators needed to account for the electroweak 
absorptive parts that arise in top pair production or annihilation. In particular, we show 
that up to NNLL order there are no contributions from interference effects originating from 
(ultrasoft) gluons carrying momenta of order m t v 2 . In Sec. lllll the electroweak absorptive 
parts of the matching conditions for the (e + e~)(t t) effective theory operators relevant for the 
bW + bW~ final state are computed. In Sec. llVl the resulting NNLL corrections for the total 
cross section are determined and the renormalization of the (e + e~)(e + e~) effective theory 
operators needed to account for the phase space divergences is discussed. In particular, we 
compute and solve the anomalous dimensions and determine their contribution to the total 
cross section. Sec. |V] contains a brief numerical analysis and in Sec. lVII we conclude. 



II. POWER COUNTING AND MATCHING CONDITIONS 

In this work we are interested in the total cross section and not in any differential in- 
formation on the top decay final states. We therefore include all effects related to the top 
quark decay as non-hermitian matching conditions of effective theory operators that describe 
the non-relativistic top and antitop dynamics and their interactions with soft and ultrasoft 
gluons. We employ gauge invariant operators, and the matching conditions are computed 
for on-shell external lines. This allows to maintain gauge invariance in a transparent way. 

To illustrate the power counting needed to classify the order at which these electroweak 
effects can contribute let us recall the matching conditions for the bilinear quark field opera- 
tors. They are obtained by matching top or antitop 2-point functions in the effective theory 
to those in the full electroweak and QCD theory. The result up to NNLL order reads 

£ ^ S) = X>J(*){^ - {J ^ L + £3 + (l - £2) - 8m t y p {x) 

+ - Xp) , (2) 

and includes the terms shown in Eq. (JTJ). Here, the fields ip p and \p destroy top and antitop 
quarks with momentum p, = (D°, — D) = <9 M + igA^ is the ultrasoft gauge covariant 
derivative and T t is the top quark width defined at the top quark pole. At order g 2 , g being 
the SU(2) gauge coupling, the width has the form 

a\V tb \ 2 m t 2 

r t = - 1 ^-(i- ar ) (1 + 2*), (3) 

where s w (c w ) is the sine (cosine) of the weak mixing angle, a the fine structure constant 
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and 



-4- W 

mf 

We use the usual v -counting D° ~ m t v 2 ~ T t ~ m t g 2 , which leads to the scaling relation 

v ~ a s ~ (7 ~ g (5) 

for the SU(2) and U(l) gauge couplings g and g' . Because the weak mixing is of order one 
we apply the same counting to the SU(2) and U(l) gauge couplings. The term 5m t is a 
residual mass term of order v 2 that arises if a threshold mass scheme jfjf is used. The LL 
order terms in Eq. (J2J) lead to the top/antitop propagator 

(6) 



P° - P 2 /(2m t ) + iT t /2 - 5m t 



Although the exchange of time-like ultrasoft A gluons contributes at LL order according to 
the f-counting, their contribution at LL order can be removed from the particle-antiparticle 
sector of the theory by a redefinition of the top and antitop fields related to static Wilson 
lines (221 . |23j]. 4 The time dilatation term cx r t (p 2 /2m 2 ) originates from the momentum- 
dependence of the full theory spinors and contributes at NNLL order. 

The ti pair is produced by an electroweak process. As long as electroweak effects are 
only treated at leading order in the f-expansion it is sufficient to describe ti production 
by a bilinear quark-antiquark current. However, as shown below it is necessary to include 
the initial-state (e + e~) fields to ensure electroweak gauge invariance at sub leading order in 
the w-expansion. The dominant operators that have to be used describing ti spin-triplet 
production have the form 

O v>p = [e 7j e]O p>1 , (?) 
Oa, p = [e 7 i7 5 e]C£i, (8) 



where 



^ p a 3 {io 2 ) X % ■ (9) 



They give the contribution AC = J2 P (CyOv,p + CaOa, p ) + h.c. to the effective theory 
Lagrangian where the hermitian conjugation is referring to the operators only. The index 
j = 1, 2, 3 is summed. The corresponding operators describing ti annihilation are obtained by 
the hermitian conjugation. Since in this work we only focus on the electroweak effects related 
to the top quark decay and, in particular, neglect QED radiative corrections (including QED 
binding and the beam effects mentioned above) the electron and positron fields act like classic 
fields in the effective theory. 5 In a more complete treatment of electroweak effects, however, 



4 In an explicit computation, quark pair production and quark-antiquark scattering diagrams involving the 
time-like gluons cancel at LL order. In Coulomb gauge all leading order diagrams with time-like gluons 



are dimensionless and individually zero in dimensional regularization ^ , . 
5 For a treatment of nonrelativistic QED effects within vNRQED see Ref. |2J] . 
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their interactions with photons have to be accounted for |20j. The leading order matching 
conditions of the operators CV/Ap obtained from the full theory Born diagrams with photon 
and Z exchange are of order g 2 and read 



a 



born / 
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[V 
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m 2 {4cl 
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(10) 
(11) 



where v is the vNRQCD renormalization scaling parameter. 6 Except for the QED beam 
effects, which we will not consider here, the electroweak effects in Eqs. (1311011 1|) are the 
only ones at leading order. In particular, at this order there are no electroweak effects 
contributing to the Coulomb potential acting between the ti pair, 



C 



pot 



p,pf w v > 



(12) 



where Vc S \u) = —AnCpC( s {m t v) is the LL Coulomb Wilson coefficient for a color singlet 
heavy quark pair. 

For the NLL order approximation it has been frequently stated that there are no new 
operators that can contribute and that power counting tells that we only need to consider 
0(a s ) QCD corrections to the LL matching conditions in Eq. @ to account for all elec- 
troweak effects (25I . 26[. However, it was also noted in {19] that the mismatch between the ti 
phase space in the full and the effective theory leads to additional NLL matching corrections 
for unstable top quarks. We come back to the role played by these specific NLL contributions 
in Sec. llVI and ignore them for the following considerations. Thus, concerning NLL effects 
related to the top decay only the one-loop QCD corrections to the on-shell top decay width 
have to be accounted for [27]. In particular, there are no QCD interference effects from 
gluon radiation off the top/antitop quark or its decay products [25, 26]. In our approach, 
which only aims at the total cross section, one can show that such QCD interference effects 
do not even contribute at NNLL order in the non-relativistic expansion. 

To discuss the QCD interference contributions we have to consider ultrasoft gluons, which 
carry momenta of order mtv 2 ~ T t , because they can interact with a resonant top quark 
without kicking it off-shell. For the time-like A gluons we already mentioned that their 
leading interaction with the quarks can be removed by a field redefinition related to static 
Wilson lines 0, H^]. Moreover, QCD gauge invariance ensures that the dominant elec- 
troweak matching corrections to the A interaction vertex vanish because we can set the 
ultrasoft gluon momentum to zero. Radiative corrections can, however, generate an anoma- 
lous interaction in analogy to the g — 2 in QED. Yet this interaction is suppressed by a factor 
1 / m 2 and cannot contribute to the cross section matrix elements at NNLL level without even 
having accounted for additional powers of the coupling constants. It remains to discuss the 



6 In vNRQCD the renormalization scales for soft and ultrasoft fluctuations, /15 and fiu, are correlated 
through the heavy quark equation of motion, = fx^/mt- The correlated running from the hard scale 
down to the soft and ultrasoft scales is described by the dimensionless scaling parameter v defined by 
/is = TritV and \x\j — m t v 2 . Thus v = 1 corresponds to the hard matching scale. 
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space-like ultrasoft A gluons which couple to the quarks with the p.A/m t coupling. In pure 
QCD space-like ultrasoft gluon exchange contributes to the cross section matrix elements at 
NNLL order and also to the renormalization group running of the ti production operators at 
the NLL level [1[ through operator mixing. Accounting for the g 2 -suppression from an addi- 
tional electroweak loop correction to the interaction vertex then also leads to a contribution 
beyond the NNLL order. 

At NNLL order let us first consider whether one needs to account for any operator in 
addition to those present already at the LL level. Ultrasoft gluon interactions have been 
discussed above. Soft gluon operators first contribute at the NLL level in pure QCD (for 
example as corrections to the Coulomb potential H^]), so g 2 electroweak corrections 
to these operators are beyond NNLL order. Such corrections can also not contribute at 
NNLL order through mixing since the Coulomb potential does not cause UV divergences. It 
remains to discuss four-quark operators. Because an electroweak loop would require a factor 
g 4 in addition to the 1/m 2 suppression for dimensional reasons, such an operator could not 
contribute at NNLL order either [30l |. 

It remains to discuss g 2 corrections to the operators contributing at the LL level. For 
the bilinear quark operators in the effective theory Lagrangian one obtains the terms shown 
in Eq. (j2J). Concerning the instability of the top quark only the time dilatation correction 
is obtained, and the 0(a 2 ) and one-loop electroweak corrections to the on-shell top quark 



width |31|, |32| have to be accounted for. The 0{a 2 ^) correction to the top width is easy to 
implement together with the Born and one-loop QCD results in the width term T t and will 
not be discussed further in this work. On the other hand, for the Coulomb potential all 
dominant g 2 corrections cancel due to SU(3) gauge invariance because in the first approx- 
imation one can neglect the gluon momentum flowing into the vertex correction |33l . liil ] . 
The mechanism is equivalent to the gauge cancellation discussed above for the time-like A 
gluon. The order g A matching conditions of the production operators Oy, P and Oa, p , on 
the other hand, do not cancel and have to be determined from matching to the one-loop 
Standard Model amplitudes for the process e + e~ — > 7, Z —>■ ti. At this order the matching 
computation can be carried out for an on-shell top-antitop pair at rest. The full set of one- 
loop electroweak corrections was determined in Ref . ■ F° r the examinations in this work 
we have to account only for the bW + and bW~ cuts because, as shown in Sec. llVl only these 
cuts are relevant for the bW + bW~ final state that can interfere with top pair production. 
Because in Ref. jl8(| only the sum of all contributions was presented, including the bb and 
W + W~ cuts, we rederive the results for the bW + and bW~ cuts in the next section. 

III. ABSORPTIVE MATCHING CONDITIONS 

The top pair production diagrams in the full theory that need to be considered to de- 
termine the absorptive bW + and bW~ cuts are shown in Fig.0 The external (on-shell) top 
quarks can be taken to be at rest. The results for the cuts in the full theory amplitude have 
the form 



A 



(A;')7 M (< W ' abS + < W ' abS 75)« e -(fc) u t {p)lMv) , (13) 
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(b) 



(c) 





FIG. 1: Full theory diagrams in Feynman gauge that have to be considered to determine the 
electroweak absorptive parts in the Wilson coefficients C\ and Cy related to the physical bW + 
and bW~ intermediate states. Only the bW + cut is drawn explicitly. 



where k + k' = 2p = (2m t , 0) and 

■,bW,abs 1 '■' lT \^tb 
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2 8si 
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\2m 2 t s 1 w x{Ac 2 w -x)(l + x) 



3a;(l + x) ( x — 4 



4 s 2 



+ g t (l-a;)(l + 2a;)(l + x + x 2 



-(1 + 12a; + 9a; 2 + 2a; 3 ) + -Ar(2 + 41a; + 28x 2 - x 3 + 2x 4 ) 
2 Ss± 



(14) 



The results for the charge conjugated process describing top pair annihilation, on the other 
hand, read 



A 



n, a^^Cy^ + tcT'^i^v^ik')] [v- t (p)^Mp)}- (is) 
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We used the cutting equations to obtain the result and checked electroweak gauge invariance 
by carrying out the computation in unitary and Feynman gauge. In both cases we computed 
the cut W lines with physical polarizations as well as with unphysical ones including also 
the charged Goldstone exchange. We note that the contributions that arise from off-shell 
corrections in the top self-energy graphs are necessary for electroweak gauge invariance. 
Since the bb and W + W~ cuts lead to different distinct phase space factors, we found that 
it is possible to identify the results also from the formulae given in 

It is an important fact that the sign of the imaginary part of the amplitude does not 
change in the charge conjugated amplitude. As for the quark field bilinear terms discussed 
before in Eq. (J2J) this is related to the unitarity of the underlying theory. It is straightforward 
to match the amplitudes for the operators CV/Ap anc ^ ®v/a p ^° ^ ne ^ un theory results in 
Eqs. f|13|) and (jl5j) . The resulting matching conditions at the hard scale read 

L>v\V = -U = W +«tyy , 

C A {y = 1) = C b A ovn + i C^ W ' abs , (16) 

where we have included also the Born level contributions from Eqs. (I10J) and pip. We 
emphasize again that these matching conditions are valid for the operators Oy/A,p and 
Oy, Ap . In a full treatment of electroweak and QCD effects the coefficients Cy/A also include 
the real parts of the full set of electroweak one-loop diag rams indicated in Figs.^ and the 
QCD matching corrections known from previous work |l6l 35l 36^ . These corrections lead to 



an energy-independent multiplicative modification of the cross section normalization which 
is, however, not subject of the investigations in this work. The results for the real parts of 
the full set of electroweak one- loop diagrams were given in [l8[ . 



IV. TIME-ORDERED PRODUCT AND RENORMALIZATION 



Using the optical theorem the NNLL order corrections to the total cross section that 
come from the absorptive one-loop electroweak matching conditions for the operators Oy/A,p 
and from the time dilatation corrections can be computed from the imaginary part of the 
(e + e~)(e + e~) forward scattering amplitude, 



cr tot ~ -Im 
s 



C' v {v) + C\{v)) L lk A lk 



where (k + k' = (y/s, 0) and e = k/|k|) 

- k + (*OV(7 5 K-(*0l k-(A07 fc (7 5 K+(AO 



Ik 



e+spins 



= ~(k + k'f (5 lk - e l e k ) 
is the spin-averaged lepton tensor and (q = (y^ — 2m t , 0)) 
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2N c S lk G°{a,v,m t ,v) 



(17) 



(18) 



(19) 
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is the time-ordered product of the tt production and annihilation operators 3 pl and 
Cp,i u 16]. In dimensional regularization the result reads 



Aa, 



r.i 

tot 



2^Im 



2% [ C£ orn C bW ' abs + C h A oin C bW ' abs ] G°{a, v, m u v) 



born" 
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2 + (C h A °™) 2 ]5G° r (a,v,m t ,u) 



(20) 



where a = —Vc S \is) / Air = CFd s {m t h'). The term G° is the zero-distance S-wave Green 
function of the non-relativistic Schrodinger equation which is obtained from the LL order 
terms in the Lagrangian shown in Eqs. (0) and (fT2*|) . In dimensional regularization it has the 
form If 



a,v,m t , v) 





In ( 
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where 



i/s - 2(m t + 5m t ) + ir t 



(21) 



(22) 



\/s being the c. m. energy. The term SG^ represents the corrections originating from the 
time dilatation correction in Eq. (J2J) and reads 



<5Gp(a, v, m tl v) 



2m t 



v d d 

1 + -tt + a— 

2 ov da 



G°(a, v, m t , v) 



(23) 



Note that the Wilson coefficients Cy/A do not run at LL order, so only the matching condi- 
tions at v — 1 appear in Eq. ()20|) . 

It is straightforward to check that the terms proportional to C V ^ A A s in Eq. (J20)) are in 
agreement with the full theory matrix elements from the interference between the double- 
resonant amplitudes for the process e + e~ — > tt — > bW + bW~ (Fig.|2t) and the single-resonant 
amplitudes describing the processes e + e~ — > t + bW~ —> bW + bW~ and e + e~ — > bW + t — > 
bW + bW~ (Figs.|21>i) in the tt threshold limit for m t — > oo. Note that diagram (a) dominates 
in the non-relativistic limit due to two resonant top/antitop lines, while diagrams (b-i) 
are v 2 -suppressed having only one resonant top/antitop line. Diagram (a) also contains 
a subleading v 2 -suppressed contribution that has to be accounted for. Diagrams with no 
top/antitop line are suppressed by v A and do not need to be considered. This also means 
that pure background diagrams containing no intermediate top quark can be neglected at 
this order. We also note that to find literal agreement between full and effective theory 
matrix elements one has to replace the ie terms in the resonant full theory top propagators 
by the Breit-Wigner term im t T t /2. The circle shown in Fig.|2k represents the QCD form 
factors for the tt vector and axial-vector currents. In the non-relativistic limit they reduce 
to the insertions of Coulomb potentials described by the higher order terms in Eq. (j21j) . Due 
to the cancellation of the QCD interference effects caused by gluons with ultrasoft momenta 
there are no further QCD corrections in the non-relativistic limit. 
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FIG. 2: Full theory Feynman diagrams describing the process e + e~ — ► bW + bW~ with one or two 
intermediate top or antitop quark propagators. The circle in diagram (a) represents the QCD form 
factors for the tt vector and axial- vector currents. 



An interesting new conceptual aspect of the corrections shown in Eq. (|20|) is that they have 
UV 1/e-divergences that arise from a logarithmic high energy behavior of the top-antitop 
effective theory phase space integration for matrix elements containing a single insertion 
of the Coulomb potential. In the forward scattering amplitude these UV divergences arise 
because the imaginary parts of the matching conditions of Eqs. lead to a dependence on 
the real part of G° (see Eq. (J21|) ). In the full theory this logarithmic behavior is regularized 
by the top quark mass. While phase space logarithms are known in the literature and can 
be resummed with renormalization group techniques j33], the divergences here are specific 
since they would not exist if the top quark were approximated as being stable. In particular, 
the UV divergences from the time dilatation corrections arise from the Breit-Wigner-type 
high energy behavior of the effective theory top propagator in Eq. (JBJ which differs from 
the one for a stable particle. Likewise, the interference effects described by the absorptive 
electroweak matching conditions for the operators Oy v and Oa, p would not have to be taken 
into account if the top quarks were stable particles. UV divergences of the same kind have 
already been observed and described before in the NNLL relativistic corrections to the S- 
wave zero-distance Green function if the unstable propag ator in Eq. ® is used [3 H H^. 
For the P-wave zero-distance Green function, which is generated by tt production through 
an axial-vector current from the Z-exchange and that contributes only at NNLL order, a 
similar UV divergence arises already at leading order in the non-relativistic expansion. Like 
for the case of the time dilatation corrections these divergences originate from the modified 
high energy behavior of the unstable top propagator. We believe it is evident that these 
divergences do not represent a deficiency of the effective theory, because the concept of 
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separating resonant and non-resonant fluctuations appears to be the only practical way to 
make systematic predictions involving unstable particles. Thus these UV divergences should 
be handled with the renormalization techniques known from effective theories for stable 
particles. The only difference is that the renormalization procedure will involve operators 
having non-hermitian Wilson coefficients. 

The operators that are renormalized by the UV divergences displayed in Eq. (|2U|) are the 
two (e + e~)(e + e~) operators 



6 A = -[e^^e] [ej^^e] 



(24) 
(25) 



which give_the additional contribution AC = CyOy + CaOa to the effective theory La- 
grangian, Cy/A being the Wilson coefficients. Because in this work we neglect QED effects, 
the electron and positron act as classic fields and therefore Cy and Ca run only through 
mixing due to UV divergences such as in Eq. ()20|) . Since only the imaginary parts of the 
coefficients CW4 can contribute to the total cross section through the optical theorem we 
neglect the real contributions in the following. Using the standard MS subtraction procedure 
the (non-hermitian) counterterms of the renormalized Oy/A operators read 



5CV 



N c m 2 t 



32yr 2 e 
+ i 



(C 



bin']] ,2 J t _|_ 2£<born£<bW,abs 



327r 2 e 
N r m 2 



born ^2 
V 



m t 



m t 



(2c 2 (u)-l^ s \u) + Vi s h 



V) 



+ ^^(^ x ) 2 -v c w ( 



48vr 2 e 



5C 



A 



N c m 2 
327r 2 e 

N c m 2 



I >uni {1 L t ^ 2^born^(bW,abs 



(C7 rn ) 



+z I -^(c h A n — \ho 

32vr 2 e v A 1 mA\ 



m t 
m t 



v) 



+ i 



N c m 2 
487r 2 e 



T) 2 — vi a) (u) 

m t 



(26) 



where the respective first term on the RHS's subtract the 1/e divergences shown in Eq. (|20|) 
and the other terms account for the UV divergences in the P-wave Green function and 
the NNLL corrections of the S-wave Green function computed in Ref. Here, C^, A 

are Born level Wilson coefficients of operators describing top pair production in a P-wave 
(originating from pure Z exchange), is the color singlet coefficient of the potential 



p2 _|_ p /2 )/(2m i (p — p') 2 ) |9| and c 2 the coefficient of the p 2 -suppressed S-wave production 
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current ^6|. The explicit formulae read 



/'"fax 



an 



C 



ax 
.4 



mj(Ac 2 w - x) 



an 



Qe + 



4s?„ 



4s 2 w m 2 {4c 2 w -x 
Vl s \v) = -AnC F a s (m t )z 
1 8C F 



l + ^ln(2-,) 



c 2 (^) 



In 



2; = 



6 3/5 V2-2 
a s (m t ) 



(27) 



The resulting renormalization group equations for the Wilson coefficients Cy/A have the 
form 



dCyju) 

din v 



^mf f (C bor„ )2 £* ( 2c2 (i/)V<'V) + V«(f)) + 2C born C bW ' abs V c (s V) 



(ilnz/ 



N c m 2 t 



8tt 2 
iV c m 2 



(^ orn ) 2 — (2c 2 {v)V^\v) + V r (s V)) + 2C born C5 w ' abs V c (s) ( 



r, 



(28) 



and the solutions for scales below m t [y < 1) read 



C v (v) = C v {l)+i 



2N c m 2 C F 
3/3 



- 1 - * , o^yborii/^(bW,abs 



ln(z) 



4C F r\ (c , orn)2 ^ + 4(c A + 2CV) r\ (C b 0m)2p(z) 



C A (z/) = C A (l)+i 



2N c m 2 C F 



m t 



3(3 



(cTt + {CaY ) — + 3C7 ra c^ 



\n{z) 



where 



A) m < 



— (C A OTn ) 2 \n 2 (z) + 



0o 



m t 



7f2 _ 1 

12 ~ 2 



ln 2 2 + ln21n(z) - Li 2 (|) 



(29) 



(30) 
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and the CV/a(1) ar e the hard matching conditions, which are presently unknown. Finally, 
the contribution of the operators Oy/A to the total cross section reads 

(31) 

Parametrically Acr^f is of order g 6 . Comparing this with the LL cross section which counts 
as g 4 v ~ g 5 we see that it constitutes a NLL contribution. This is also evident from the 
fact that the corresponding UV divergences were generated in NNLL order effective theory 
matrix elements. The correction A<j to ' t is energy-independent, but it is scale-dependent 
and compensates the logarithmic scale-dependence in the NNLL order matrix elements. As 
mentioned before, the matching conditions Cv/a{ v — 1) ar e presently unknown and we 
therefore set them to zero in the numerical analysis presented below. We note, however, 
that it was shown in ^| that the difference between the full theory phase space (which is 
cut off by the large, but finite mt) and the effective theory phase space (which is infinite in 
the computation of the forward scattering amplitude) contributes to Cv/a{v = 1) and also 
represents a NLL effect. 



Atr. 



r,2 

tot 



Im 



C v + Cj 



V. NUMERICAL ANALYSIS 

In Fig. El we have plotted Aa^t an d A(T t r ' t 2 in picobarn in the IS mass scheme [l^, |3^| 
for Mis = 175 GeV, a = 1/125.7, s 2 w = 0.23120, V tb = 1 and M w = 80.425 GeV with the 
renormalization scaling parameter v = 0.1 (solid curves), 0.2 (dashed curves) and 0.3 (dotted 
curves). The divergences in Aa^'t are subtracted minimally. For the QCD coupling we used 
a s (Mz) = 0.118 as an input and employed 4-loop renormalization group running. Note that 
in the IS scheme 8m t = Mis(vj (^)/4tt) 2 /8. In Fig.EK the sum of Acr^t an d Acr^ 2 is shown 
while in Fig.Eb both contributions are presented separately. For the top quark width we 
adopted the value T t = 1.43 GeV. 7 We find that the sum of the corrections is negative and 
shows a moderate v dependence. Compared to the most recent NNLL QCD predictions for 
the total cross section [17| the corrections are around —10% for energies below the peak, 
between —2% and —4% close to the peak and about —2% above the peak. Their magnitude 
is comparable to the NNLL QCD corrections. Interestingly, they partly compensate the 
sizeable positive QCD corrections found in [l3|, [I?} . The peculiar energy dependence of the 
corrections, caused by the dependence on the real part of the Green function G°, also leads 
to a slight displacement of the peak position. Relative to the peak position of the LL cross 
section one obtains a shift of (30, 35, 47) MeV for v = (0.1, 0.2, 0.3). This shift is comparable 
to the_expected experimental uncertainties of the top mass measurements from the threshold 
scan 



VI. CONCLUSION 

We have determined electroweak corrections to the NNLL top pair threshold cross section 
in e + e~ annihilation related to the instability of the top quark. Our approach is closely 



7 In a complete analysis of electroweak effects the top quark width depends on the input parameters given 
above and is not an independent parameter. For the purpose of the numerical analysis in this work, 
however, our treatment is sufficient. 
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FIG. 3: The corrections Ao&jf and A<7 t r o ' t 2 in p6 for M 1S = 175 GeV, a = 1/125.7, 4 = 0.23120, 
V tb = 1, M w = 80.425 GeV, T t = 1.43 GeV and v = 0.1 (solid curves), 0.2 (dashed curves) and 0.3 
(dotted curves) in the energy range 346 GeV < y/s < 354 GeV. Panel (a) shows the sum of both 
corrections and panel (b) the individual size of Ao" to ' t 1 (energy-dependent lines) and A(7 t0 f (straight 
lines) . 

related to the treatment of absorptive processes in the optical theory. It includes the effects 
of the top instability by accounting for the absorptive parts in electroweak corrections to the 
effective theory matching conditions that are related to the observable bW + bW~ final state. 
The matching conditions render the NRQCD Lagrangian non-hermit ian, but they allow 
for the determination of the total cross section using the e + e~ — > e + e~ forward scattering 
matrix element and the optical theorem. We have shown that the absorptive parts of the 
electroweak matching conditions for the tt production and annihilation operators describe 
the interference of the double-resonant amplitude for e + e~ — > tt — > bW + bW~ with the 
single-resonant (and v 2 -suppressed) amplitudes for e + e~ — > bW + t — > bW + bW~ and e + e~ — > 
tbW~ — > bW + bW~ . We have also shown that at NNLL order there are no further interference 
effects caused by the exchange and radiation of ultrasoft gluons. 

The novel feature of the NNLL corrections is that they lead to new UV divergences. 
These divergences originate from the logarithmic high energy behavior of the tt phase space 
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in the effective theory forward scattering amplitudes, that is caused by the interferences 
and by the modified propagators of an unstable top quark. The divergences renormalize 
(e + e~)(e + e~) operators that contribute to the forward scattering amplitude already at NLL 
order. The corrections determined in this work slightly modify the cross section shape and 
are comparable to the known NNLL QCD corrections. The size of the corrections shows 
that a complete treatment of all NNLL electroweak effects is desirable. 
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